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m . 

f*^ , We derive the electrical and thermal conductivities of the d-density-wave state in an external 

■ magnetic field B in the low-temperature regime and in the presence of impurities. We show that in 

' the zero-temperature limit, T — > 0, the Wiedemann-Franz (WF) law remains intact irrespectively of 

the value of the applied field and chemical potential \i as soon as there is scattering from impurities. 
For finite T < the WF law violation is possible and it is enhanced by the external field. Depending 
on the relative values of B and T the electrical conduction may dominate the heat conduction and 
vice versa. For fi^> T the WF is restored even in the presence of the external field. 
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I. INTRODUCTION 



A recent experiment of Hill et al. [| that measured electrical and thermal conductivities of the optimally doped cop- 
per oxide superconductor Pr 2 _: E Ce :E Cu04 (PCCO) in its normal state found striking deviations from the Wiedemann- 
Franz (WF) law. The electron-doped PCCO compound with T c ~ 20K near the optimal x — 0.15 doping is the analog 
of the hole doped material La2- a: Sr : j;Cu04_j / (LSCO), which is also optimally doped at approximately x — 0.15 with 
a maximum T c ~ 40K. Lower T c in PCCO makes it particularly convenient from a technical point of view since an 
external magnetic field is necessary to destroy superconductivity to study its low temperature transport properties in 
I . 1 the normal state. While for LSCO the upper critical field is ~ 50T, in the optimally doped PCCO the superconduc- 
tivity is already destroyed at 8T pj. The choice of overdoped compounds that can be driven into the normal state 
by a relatively small magnetic field ~ 13T is wider and the hole overdoped system Tl2Ba2Cu06+5 (Tl-2201) with 
T c = 15K was studied very recently by Proust et al. 0. They verified that in the overdoped Tl-2201 WF law holds 
perfectly. 

The WF law is one of the basic properties of a Fermi liquid, reflecting the fact that the ability of a quasiparticle 
to transport heat is the same as its ability to transport charge, provided it cannot lose energy through collisions. 
t— ( ■ The WF law states that the ratio of the heat conductivity k to the electrical conductivity a of a metal is a universal 
CN | constant: 

, where ks is the Boltzmann constant, e is electron charge, and Lq = 2.45 x 10 8 W£!K 2 is Sommerfeld's value for 
the Lorenz ratio L = n/(aT). To be more precise, one should also specify the temperature range where the WF law 
holds. Strictly speaking this law is proven only in the limit T — > and for a small concentration of impurities [3]. In 
a less strict sense it is often implied that the WF law is valid for fcgT -C /i, where /i is chemical potential. In this case 
one can provide some qualitative arguments Q that if the scattering from impurities does not strongly depend on the 
energy the WF will still be valid. Moreover, for fcgT <C /i and elastic electron scattering these arguments can also 
be extended to the case of a strong (Hujl 3> fc^T, tkvL S> T) quantizing magnetic field [5j, where hu>L is the distance 
between Landau levels and T is the impurity scattering rate. Thus, if these conditions are fulfilled, any violation of 
the WF law would indicate that there is a breakdown of the picture based on the Fermi liquid theory. 

One of the possible theoretical interpretations of the WF law breakdown is that the quasiparticle fractionalizcs 
into separate spin and charge. This separation can be investigated using various models and approaches, so here we 
mention only the most recent studies done in the context of the WF law violation in cuprates. Specifically, the WF law 
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violation has been studied for a large- N limit of the t — J model on the square lattice Jof . Another examination of the 
WF law was done by Yang and Nayak (YN) and also by Kim and Carbotte (KC) jg] within the phenomenological 
d-density-wave (DDW) picture. 

The DDW scenario proposed in Ref. 9] is based on the assumption that the pseudogap phenomenon ^(j m high-T c 
cuprates is the result of the development of another order parameter called the DDW order [Tj that has d-wave 
symmetry: 

( C t(t,k+Q)cv(*,k)) =i$ Q /(k)<W (1.2) 

where s, s' = ± is the spin index, /(k) = cos k x a— cos k y a, Q = (7r/a, 7r/a), and a is the lattice constant. In comparison 
to the d-wave superconducting order the DDW order parameter appears to be rather exotic because it breaks the 
time reversal, translation by one lattice spacing and rotation by n/2 symmetries, but respects any combination of two 
of these. The underlying reason for this breaking is that there are countercirculating currents in the DDW ground 
state. The schematic phase diagram for the DDW pseudogap scenario is shown in Fig. ^ For a hole doping larger 
than the critical value x c w 0.19 [l^, the DDW order is presumed to disappear, whereas for the underdoped system 
DDW order coexists with d-wave superconductivity (dSG). 

One of the unusual features of the DDW state is that for a half-filled system the chemical potential of the nodal 
quasiparticles participating in the electrical and thermal transport can be small or even zero, i.e., < ksT, which 
violates the usual conditions of the WF law validity. Indeed, exactly in the limit \i = the WF law is strongly violated 
in the extremely clean limit t 7]. There is no WF law violation in the T = limit for finite \i and/or r |3,|8j. For finite 
temperatures the WF violation depends on the impurity scattering 8]: in the Born limit (for a constant impurity 
scattering rate) there is no change in the WF law, but in the unitary limit for the frequency dependent scattering 
rate the WF law is violated, but only for smaller that the DDW gap. When \i is increased sufficiently, the Lorenz 
number becomes approximately equal to its conventional value and its temperature dependence is small. 

While in general the validity of the DDW pseudogap scenario is still questionable, it is important to scrutinize all 
its theoretical consequences. One of the opportunities is to study possible WF law violations using the DDW model, 
so that these results can be compared with the experimental results of Refs. 0. Although this investigation to 
a large extent had already been done by YN and KC 0, both these papers do not consider the presence of the 
external magnetic field which is an essential ingredient of the experiments 0, • Indeed, in these experiments the 
DDW state appears in the magnetic field in the underdoped regime at low T when the superconductivity is destroyed. 

In this paper we study the WF law for the DDW model including the external magnetic field for a constant 
impurity scattering rate, paying special attention to the regime < ksT where the violation of the WF law is 
expected. Although there are some formal similarities between the DDW and d-wave superconducting states, the 
present problem is much simpler because there are no vortices and no superflows surrounding them in the DDW 
state. The U(l) symmetry in the DDW state also remains unbroken, so that the electromagnetic field enters into the 
effective low-ener gy t heory in the same way as in QED2+1. Thus we avoid the famous problem of Dirac Landau-level 
mixing (see Refs. |l3tll4j'). 

We begin by presenting in Sec.[n]the model DDW Hamiltonian and rewriting it using two-component spinors. Since 
there is a discrepancy between the expressions for the electrical current used by YN and KC |g, we also derive 
the electrical current operator and show that in the nodal approximation it reduces to the current of Ref. . Then 
we obtain the low-energy effective action (Lagrangian) for the nodal quasiparticles in the presence of an external 
electromagnetic field, consider its discrete symmetries, and rewrite the Lagrangian in the Dirac form. In Sec. IIIII 
we consider the Green's function of nodal quasiparticles in an applied field (the calculational details are given in 
Appendix [SJ . In Appendix [5] using this Green's function, we derive a generalized polarization bubble that can be 
applied to the calculation of electrical and thermal conductivities. In Sees. II Vl and Ivl we obtain general expressions 
for electrical and thermal conductivities in the external field. These general expressions are further analyzed in detail 
by analytical and numerical methods. In Sec. I VII we investigate the implications for the WF law that follow from the 
results presented in the previous sections. Conclusions are discussed in Sec. IVIII where we give a concise summary of 
the obtained results and compare them with experiment. Some useful integrals are given in Appendix IO 

II. d-DENSITY WAVE HAMILTONIAN AND ITS EFFECTIVE LOW ENERGY FORM 

We start with the mean-field Hamiltonian for the DDW state 0, El 113 

# DDW = / ^Mk)-^cl(k)c s (k) +l D(k)ct(k)c s (k + Q)}, (2.1) 
Jbz 

where the single particle energy is e(k) = — 2t(cosk x a + cos k y a) — At' cos k x a cos k y a with t, t' being the hopping 
parameters, [i is the chemical potential, D(k) — ^(cos k x a — cosk y a) is the ci-density-wave gap and Q = (n/a,n/a) 
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is the wave vector at which the density-wave ordering takes place. The integral is over the full Brillouin zone. 
Throughout the paper h — ks = c = 1 units are chosen, unless stated explicitly otherwise. 

Due to the presence of the factor i in the second term of the Hamiltonian 1(2.1)) . it breaks time-reversal symmetry, 
where the time-reversal operation is defined as |16| 

C(k) -» Tcs^T- 1 = isc_ s (-k), 4(k) -» TctO^T- 1 = -iscl s (-k) (2.2) 

and the time-reversal operator T is unitary and antilinear. Nevertheless the Hamiltonian 12.1fl preserves a combined 
symmetry of the time reversal and a translation by one lattice spacing T = 7~T a : 

c s (k) -» Tcs^f- 1 = te-^c^(-k). (2.3) 

Since the order parameter also breaks translational symmetry by one lattice spacing, it is convenient to halve the 
Brillouin zone as shown in Fig. [21 and form a two-component electron operator 

Xs(t,k)= ( Cs ^+ Q) ). xU*,k) = ( C t(t,k) ct(t,k + Q)). (2.4) 



Then the mean-field Hamiltonian (|2.1|l in terms of x becomes 

Xs(t,k), (2.5) 



d 2 k 



Vrbz (2tt) 2 



i[e(k) + e(k + Q)] - » + ± [e(k) - e(k + Q)]<7 3 - D(k)a 2 



where Ui are Pauli matrices and the integral is over the reduced Brillouin zone (RBZ). The symmetry operation (|2.3|l . 
for example, is 

Xs (t,k) -> fxsit^f- 1 = -e~ tka a 2 x~s{-t, k - Q). (2.6) 

In what follows we will consider the more simple case t' = 0, so that after employing the nesting property e(k + Q) = 
— e(k) the Hamiltonian (|2.5|l takes a simple form 

HDDW = ! T^X^k) [H (k) - fi\ Xs (t,k), (2.7) 

where 

i?o(k) =£(k)a 3 -^(k)a 2 . (2.8) 



A. Electrical and energy current operators 

As we already mentioned the expressions for the electrical current used in Refs. and are different. Although 
quantitatively this difference for high-temperature superconductors is expected to be rather small [see after Eq. Q2.21[l ]. 
the extra term in the current present in Ref. |7j is of qualitative importance, so we need to discuss its origin. 

To find out electrical current operator j for Hamiltonian 1)2. 1|) we apply the method similar to that of Refs. 01 an d 
where it was used to determine the heat current in the <i-wave superconducting and DDW states, respectively. 
The idea of the derivation is to make the expression for the current operator compatible with the charge conservation 
law, 



d P {t,v) 
dt 

and with the equations of motion for field operators, 

,&4(k) 



V-j(t,r) = 0, (2.9) 



dt 

Writing Eq. 1)2.9(1 in the momentum space, one obtains 

dp(t,q) 



[ C t(k), H] = -c+(k)[e(k) - M ] - 4(k - Q)zD(k - Q). (2.10) 



t)l , *q-j(*,q) = 0, (2.11) 
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where the Fourier transform of the charge density is 



For its derivative we obtain 



Wq) 

at 



BZ 



d 2 k 

(2tt)= 



pit, q) 



c|(k)c s (k + q)q 



BZ 



d 2 k 
{2-Kf 



de{k) 
dk 



4(k)c s (k + q). 



f(k)c fl (k + q+Q)iq 



8D{k) 



so that the electrical current is 

j(*,0) 



d 2 k 
d 2 k 



c t(k) Cs (k 



0e(k) 



5k 

xi(k)V(k) X .(k), 



•4(k)c s (k + Q)i 



Ok 



dD{k) 
dk 



(2.12) 



(2.13) 



(2.14) 



/RBZ ( 27r )" 

where we introduced the generalized velocity operator that depends of the Fermi, yp , and gap, V£>, velocities 



V(k) = v F (k)cr 3 + y D (k)a 2 , Vj?(k) 



9k ' 



vr>(k) = - 



ftD(k) 
9k ' 



(2.15) 



In the last equality in Eq. I|2.14|l we used the two-component form and relied on the nesting property e(k+Q) = — e(k) 
and D(k + Q) = -D(k). 

This form of the electrical current operator reduces to the current used by YN [?| after the nodal approximation is 
made. The last term of Eq. I|2.14|l . ~ vp, is due to the fact that the DDW gap in the mean-field Hamiltonian l|2.1() 
is k-dependent. The expression for the electrical current used by KC contains only one term ~ vjr. In this case, 
calculating electrical conductivity in the bare bubble approximation is not consistent with charge-current conservation, 
since such a bare vertex does not satisfy the Ward identity when quasiparticle self-energy has a nontrivial momentum 
dependence. Instead one should use a vertex that complies with the Ward identity. 

The derivation of the energy current operator from the Hamiltonian (|2.1I) is quite similar to the derivation of 
the electrical current operator (|2.14l) . so here we only outline the main steps. This current is calculated from the 
corresponding continuity equation for the energy density and the equations of motion (|2.10() (we set there /i = 
because the energy density is considered) that give 



or in the spinor form 



'RBZ ( 27r ) 

Then using again the equations of motion we finally get 



f(t, o)= [ [g(k)v F (k) + I>(k)v D (k)] ct(k)c(k), 

jbz 

= \l 7^X t (* ) k){^o(k) J V(k)}x(t,k). 



X T (w,k)V(k)x(o; + n,k). 



/RBZ (2^) 2 

One can easily see that Eq. I|2.18|l agrees with the corresponding energy current obtained in Rcf. 



(2.16) 



(2.17) 



(2.18) 



B. Effective low-energy nodal action 



The action S corresponding to the Hamiltonian (|2.7I) is 



S 



dr 



d 2 r X i(T,r)d TXs (T,r) + H DDW (r) 



(2.19) 



d 2 k 



I dT I i^xl(T,k)[d T -v + e(k)<j 3 -D(k)a 2 ] Xs (T,k), 

J jrbz { 27r r 
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where r is the imaginary time. Hence the Green's function of the DDW state reads 

k) = (ia, + ^)2_ g 2( k ) ' s ( k ) = V £-(k) + £ 2 (k), (2.20) 

where iu) — i(2n + 1)ttT is fermionic (odd) Matsubara frequency. 

Linearizing the spectrum about the four nodes N = (±7r/2a, ±7r/2a) at half-filling (/x = 0), one obtains the Green's 
function (GF) for a given node (we choose the local nodal coordinate systems as shown in Fig. |2J) 

SnoAk) = iiU + ^tV(T Va2 > W) = y/*kZ + v W, (2.21) 

where the Fermi velocity for half-filling vp = |9e(k)/9k|k=N| = 2y/2ta and the DDW gap velocity vp, — 
|9£)(k)/9k|k=N| = "^Tf^oa- It is important to stress that this way of linearization is different from the nodal ap- 
proximation used to describe dSC when the expansion is done around the nodal ("Dirac") points defined as an 
intersection of the Fermi surface e(k) — /x = and the nodal lines D(k) = 0. For the DDW case one always expands 
around half-filling e(k) = and D(k) = 0, so that for nonzero /x the nodes transform into small pockets (see Fig. 0) 
which are cross sections of the "Dirac cone" . Thus only at exactly half-filling is there a correspondence between DDW 
and dSC cases. There are the following estimates 0] of the model parameters, t = 300 meV and Dq — 50meV, which 
give the ratio vf/vd — 24. This shows that the second term of the electrical current (|2.14|) . which is proportional to 
vp, is indeed small. The value of the chemical potential corresponding to a hole doping of 10% is 01 A* = 0.91i. 
The linearized Dirac action that corresponds to the Green's function (|2.21|) can be written as 

P 

Snodc = ~ J dr J rf2r ^«( T ' r )[^(^ ~ A 1 ) - v F a 3 id x - v D a 2 idy\xs{T,v). (2.22) 
o 

Dealing with Eqs. (|2.21l) and (|2.22l) . one should not forget that all physical quantities involve the summation over the 
4 nodes present in the original Green's function l|2.20l) . In what follows we will use the nodal action (|2.22() instead 
of the original DDW Hamiltonian (|2.1|) . This approximation is well justified for T <§C Dq. As we will see below, this 
approximation and the fact that the electromagnetic field enters the theory through the minimal coupling prescription 
will allow us to apply the results of QED2+1. 



C. Discrete symmetries of the nodal Lagrangian 

Since the action (|2.22|l does not directly correspond to the common QED2+1 form used in the literature ( see 
also Ref. 20]), we present below the transformations for the discrete symmetries. Let us consider only one node and 
suppress the spin index s, so that in real time t = —ir the action S reads 



dt / d 2 rC, 



(2.23) 



with 



C 



x\t, r) ii(d t - ieA (t, r))+fi + iv F a 3 (d x 



i^Ai(t,r)) +iv D <T 2 (fy-i^4 a (t,r))] x(i,r), (2.24) 



where A u (t, r), v = 0, 1,2 is the electromagnetic field vector potential. This is a highly nontrivial fact in that the 
vector potential can be inserted in Eq. (|2.24l) using the minimal coupling prescription of QED. Only due to this fact 
can one apply the results of QED2+1 for the description of the DDW state. Moreover, this is only correct for the 
DDW state but not for the dSC state, where the last term of the Lagrangian (|2.24|) does not couple with the 
electromagnetic field. It has to be admitted that the way of inserting the electromagnetic field is not just related to 
the nodal approximation and can be traced back to the expression (|2.14|) for the electrical current that can only be 
obtained when the vector potential enters both e(k) and D(k) terms in Eq. I|2.1() . Finally we note that the symmetries 
of the Lagrangian have to be analyzed in real time where one can distinguish time and spatial coordinates. 
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1 . Parity 

In 2 + 1 dimensions, parity corresponds to inverting one axis (since the inversion of both axes would be rotation by 
7r): Py — r', where r = (x, y) and r' = (x, —y). The corresponding operation on the two-component spinor and on 
the gauge field is 

X(r, r) -» V X (r, v)^- 1 = a 3 x(r, r'), x\t, r) -> V X ^ (*, r)7> _x = X + (*, r>3, 

4>(«, r) -> PA (t, r)^- 1 = A (i, r'), Ai(i, r) -> r)^" 1 = Ai(t, r'), A a (i, r) -> 7>A 2 (i, r)^" 1 = -A a (t, r'). 

(2.25) 

One can verify that the Lagrangian (|2.24|) is invariant under this parity transformation while possible mass terms 
with <r\ and 172 break parity: 

xHt,v)a lX {t,v) -> Px t (* > r)a-ix(«,r)p- 1 = - X f (t, r> lX (i, r'); 

. , (2. 2o) 

X f (t, r)o-ax(t, r) -> V X \t, r)a 2X (t, r^ 1 = - X Ht, r')a 2X {t, r'). 



2. Charge conjugation 

Charge conjugation, 

X (t,r) ^C X {t,r)C- x = a 3 [ X Ht,r)] T , x f (*, r ) _ C ^{t, r)^ 1 - [ X (i,r)] T a 3 , 

A„(t,r) ->£M„(i,r)C- 1 - -A„(t,r), ^ = 0,1,2 

leaves the equations of motion invariant. It is easy to check that the Lagrangian (|2.24l) is invariant under this 
transformation only at half-filling (fi = 0) while away from half-filling the term with the chemical potential breaks C 

^X ] {t,v)i X {t,v) -> fj£xHt,r)i X (t,r)C- 1 = -/i X t (t,r)/x(t,r). (2.28) 
One can also check that the mass term with 02 also breaks C, whereas <j\ mass term does not break C . 



3. Time-reversal symmetry 

The operation of time reversal for the electromagnetic field is defined as 

A (t, r) -> TA (t, v)T~ x = A (-t, r), A h2 (t, r) -> TA h2 (t, r)^ 1 = -A 1;2 (-t, r). (2.29) 
We determine a matrix T responsible for the time-reversal transformation of spinors 

X (t, r) - T X (t, r)T _1 = T X (-t, r), x f (i, r) - T X f (t, r)^ 1 = x f (-i, r)T (2.30) 
demanding that the current components transform as follows 

X f (t, r)a„x(* , r) -> x Tf (-t, r^X 1 (-i, r) = X f (-*, v)Ta* v T X {-t, r) = X f (-t, r)a„x(-i, r), (2.31) 

where 

a v = (/, 0-3,0-2), a„ = (ao,-ai,-a2)- (2.32) 

The time-reversal operation (I2.3U|) can be written down using T = <r 2 and one can see that the Lagrangian C is 
invariant with respect to this symmetry 

£[ X f (i, r), X (t, r), i4 v (i, r)] -> £[ X T t (-i, r), X T (-t, r), A^(-t, r)] - £[ X f (t, r), X (t, r), A„(t, r)], (2.33) 

while the Ox and 02 mass terms would break it 

XHt, r)a h 2X(t, r) - T x f (i, r)o 1)2 x(i, OT" 1 = -x f H, r)ai, 2 x(-t, r). (2.34) 

One can check that if we use notations corresponding to those of Ref. [20|. P and T operations introduced here 
agree with the transformations given in Ref. 20]. The time-reversal symmetry operation H2.30|l for the continuum 
nodal Lagrangian should be identified with the combined symmetry operation ({2.6(1 that leaves the Hamiltonian 
(|2.5|) invariant. We note that, in principle, the time-reversal transformation should also flip the spin s — ► — s, but 
since all considered expressions assume summation over two spin components, we do not include this spin flip in the 
transformation of Dirac spinors. 

Parity and time-reversal symmetries can still be broken by the mass (gap) term ~ <j\. In particular, this mass term 
may be generated in an external magnetic field due to the so called magnetic catalysis phenomenon [2l|. 
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D. Final form of the nodal Lagrangian 

We choose the vector potential for the external magnetic field B in the symmetric Poincare gauge A = (— -f xi), 
so that the field B is perpendicular to CuC>2 planes. To consider also a possibility of magnetic catalysis [21j we add 
to the Lagrangian i|2.24|l the interaction term originating from particle-hole attraction 

Ant = \{x ] {x)a lX {x)f, x = (t,r). (2.35) 

The simplest way to treat the interaction (|2.35|) is to introduce Hubbard-Stratonovich field f{x) — — g\ (x)&ix( x )i 
so that the Lagrangian (|2.24() becomes 

C = X ^{x)[ia v D v - u W {x)\x{x) - ^M, (2.36) 

where the covariant derivatives are 

fhdt — ieAo(x), v = 0, 

v F {hd x - if A x {x)) , u=l, (2.37) 
v A (Hd v - ilA 2 {x)) , u = 2, 

and a matrices were defined in Eq. I|2.32[) . 

Finally, to simplify further calculations and to make a direct link with QED2+1 it is convenient to introduce Dirac 
conjugated spinor x(x) = x* { x )&i an d rewrite the Lagrangian l|2.36|) as 

C= x (x)[i^D v -^x)] X (x)-^-, (2.38) 

where the 7 matrices are 

7" = (<7i,-W2,ur 3 ), {7",f}=2V", ^ = diag(l,-l,-l). (2.39) 

Note that the Zeeman term, if necessary, can be added explicitly both to the original Hamiltonian (|2.1I) and the 
Lagrangian H2.38J) . Here, however, this term is neglected. 

III. GREEN'S FUNCTION IN AN EXTERNAL MAGNETIC FIELD 

To derive the transport coefficients we need an explicit representation for the fermionic Green's function in an 
external field. The calculation, following the Schwinger (proper time) approach [22|, is sketched in Appendix lAl and 
the result is (we set for convenience vf — vd — 1 and restore them when it is needed according to the prescription 
given in Appendix [5J| 

S(x - y) = exp j A?d* j S(x - y), (3.1) 
where the translation invariant part S in the Matsubara frequency-momentum representation is given by 



S(iuj, p) = — J ds exp 



,2 .,-> , „2 ta nhjegs) \ 
eBs J 



-s I A - {iujf + V l 



/J (3.2) 

x (iuj -pi-f 1 -P27 2 + A + -(P27 1 -Pi7 2 ) tanh(eBs)j ^1 + --f 1 ^ tanh(e_Bs)^ . 

Here we replaced the Hubbard-Stratonovich field ip(x) by its mean-field value A = (ip(x)}. The chemical potential /1 
has to be taken into account via the shift iu; — > iui + \i [see Eq. I|B3|I in Appendix IB]. 

As already mentioned, the gap A can be generated by an external magnetic field [2 1| . The value of the gap has 
to be determined from the minimum condition of the corresponding effective potential, see e.g. |2ll l23l l2^j . Here, 
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however, we are mainly interested in the case of A = and will concentrate on \x and T dependences of the Lorenz 
ratio L. 

The propagator l|3.2jl can be decomposed over the Landau level poles 0,113 (see also Appendix El for the details): 



S(iu;,p) = e xp (-JU5>l)« 



S n {B,iuj,p) 



where 



S n (B,iu,p) = 2(A + iu 1 °) 



P_L r 



eB 



+ -^n-l 



(iu) 2 - A 2 - 2\eB\n 



(3.3) 



(3.4) 



with P± = [1 ±sgn(ei?)j7 1 7 2 ]/2 = (1=fcti)/2 being projectors and L n , L n Laguerre's polynomials {L 1 _ l = 0). In what 
follows for convenience we take eB > 0. The spectral function A]j(ui, p) associated with the translation- invariant part 
S of the Green's function is defined as 



A D (u,p) 



1 

27ri 



S A (co - i0,p) - S r (uj + i0,p) 



(3.5) 



where the retarded, S R , and advanced, S A , Greens' functions are S R (u>+iO, p) = S(iu> — » u>+ 

S(ioj — ► uj — i0,p). Using the definition 1)3. 5f) one obtains the spectral function associated with the "Dirac" fermion 
Green's function (xx) 113. 3|l for a clean system 



A D ( W) p) = e X p(-^)^(-l 



(7°M n + A)/ 1 (p) + / 2 (p) 
2M„ 



S(u - M n ) + 



( 7 °M n -A)/ 1 (p)-/ 2 (p) 



2M„ 



where M n = VA 2 + 2eBn and 



/i(p) 



P - L " 1 2 fB 



P + L n . x ( 2^ 



/ 2 (P)= 4(pi7 X +P27 2 )ii-i (2^ 



5(w + M„) , 
(3.6) 

(3.7) 



To consider transport phenomena for a more realistic case one should introduce into the theory the effect of 
scattering on impurities. In general this can be done by considering dressed fermion propagators that include the self- 
energy due to the scattering on impurities. This self-energy in turn has to be found self-consistently by solving 
the corresponding Schwinger-Dyson equation. This scheme was in fact used in the Ref. 0, and the WF violation 
obtained by KC at finite temperatures in the unitary limit is essentially due to the nontrivial frequency dependence 
of the scattering rate T(u>) = — ImS fl (w). As mentioned in Introduction, here we choose the case of constant width 
r = T{u> = 0) = l/(2r), r being a mean free time of quasiparticles, so that the <5-like quasiparticle peaks in Eq. I|3.6() 
acquire a Lorentzian shape: 



6(w ± M n 



7T (w ± M n y + r 2 ' 



(3.8) 



Such a broadening of Landau levels was found to be a rather good approximation valid in not very strong magnetic 
fields UJ. 

Finally, we establish a link between the Green's function l|3.1|l considered here and the Green's function (|2.2U|) 
used in the preceding section. Taking into account that the average (xx } is related to the "Dirac's average" via 
(XX ) = (xx)7° one obtains for B = 0: 



G(zc J ,p) = 5(z^ m ,p;S = 0)7° = 
so that restoring vp, vd, replacing iui. 



{iuj) 2 - p 2 - A 2 ' 
H, and setting A = we recover the Green's function Q2.21|) . 



(3.9) 
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IV. ELECTRICAL CONDUCTIVITY 
A. General expression for electrical conductivity 

The frequency dependent longitudinal electrical conductivity can be calculated from the current-current correlation 
function 



IT^(iQ) =- J dre^iTrjKr, 0)^(0,0)) (4.1) 
o 

by means of the Kubo formula 

am = JmnOO(n + iO) ) (42) 

where LT^ C (ft + iO) = U cc (iQ — > ft + iO) is the longitudinal polarization [see Eq. 1B14|| ]. 

The expression for the electrical current operator for the Hamiltonian (|2.1|) was already derived in Eq. 1)2. 14[) . 
Having the current, we are in a position to calculate the current-current correlation function (|4.1() . which is given by 
the bubble diagram. The influence of the impurity vertex corrections on the transport properties of another (dSC) 
nodal system was considered in Ref. 17']. The calculation of the bare bubble polarization function is in fact similar for 
all transport coefficients in dSC and DDW nodal liquids, so that it is convenient to define a generalized polarization 
function H" that depends on the coupling parameters g, g' and the generalized velocity V(k). The polarization 
function is calculated in Appendix iBl Using the general result (|B11|) for the case of interest, g = g' = e in the limit 
il — > 0, we find that the dc conductivity of isotropic system is given by 

oo 

° = ™ 2 I 7^ [ M-n' F (Lu~ri)tr[A(u,k)V a (k)A(LU,k)V a (k)], (4.3) 

Jrbz ( 27r ) J 

— oo 

where — n' F {uo — fj,) = (1/4T) cosh~ 2 [(c<j — /i)/2T] is the derivative of the Fermi distribution and summation over the 
dummy variable is implied. 

At this point it is instructive to compare Eq. I|4.3[) with the expression for the conductivity obtained by Langer 
[Eq. (4.8) of Ref. [28j] in the bubble approximation. Langer's expression has the same form as Eq. (I4.3|l except the 
matrix- valued eV a (k) is replaced by the vertex with coinciding fermion momenta and energies 

A i (k,u)=e^-e^(k,u), (4.4) 
m oki 

where S(k, uj) is the usual self-energy [he considered the quadratic dispersion law, so that de{k)/dki — fcj/m]. In fact 
Eq. 1)4. 4f> is a direct consequence of the Ward identity. Thus knowing the self-energy one immediately obtains the dc 
conductivity in the bubble approximation. Furthermore, the bubble conductivity of Ref. [28| becomes exact in the 
zero-temperature limit. 

Since in the mean-field approximation S(k, lu) is replaced by the DDW gap 02-D(k), Langer's expression reduces to 
Eq. (|4.3|) . It is also evident from Eq. (|4.4|) that for the momentum-independent self-energy the bubble conductivity 
coincides with the conductivity calculated in the bare bubble approximation, but clearly this is not the case of the 
DDW gap. Expression l|4.3|l can also be derived in the lowest order approximation from a g eneral expression for the 
dc conductivity, with the vertex corrections taken into account, obtained by Eliashberg |29j- 

Making the nodal approximation (see Eq. I|B16|) ) we get 



2ne 2 f d 2 p 

v F v D J (2tt) 2 J r "" (4.5) 

— OC 

x (^tr [A D (w,p) 7 1 A I) (c i ;,p)7 1 ] +v 2 D tx[A D {uj, V ) 1 2 A D {uj^) 1 2 ]) 



Finally we note that one could also calculate conductivity directly working with the nodal Lagrangian (|2.38|) and 
expressing the bubble in terms of the corresponding electrical current operator 

jx = evpXs^Xs, j y = ev D Xs1 2 Xs, (4.6) 

which is also valid in the presence of an external field. This way of calculation is exactly the same as in Ref. and 
final result agrees with l|4.5|l . 
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B. Calculation of conductivity 



Straightforward calculation of the trace in Eq. 14. 51) . with Ar){ui,\i) from Eqs. (|3.6|l and l|3.7|) but with the S 
functions in Eq. (|3.6|> replaced by the Lorentzians l|3.8[l . gives 

4p2 / 2p 2 
tr[A D ( W ,p)7 1 ' 2 A D (a;,p)7 1 ^] = — cxp — 

f> f nn+m+1 [(^ 2 + M 2 + L 2 )(c 2 + Ml + L 2 ) - 4c 2 A 2 ][£»£ m _ 1 + L„_iL m ] T 32(p 2 - gfjL^L^ ( 4J ) 
^ 1 j [(w 2 + M 2 +r 2 ) 2 -4o; 2 M 2 ][(a; 2 + M^ + r 2 ) 2 -4w 2 Af 2 n ] 



n,m— 



where all Laguerre's polynomials depend on 2 ^-, and the minus sign corresponds to 7 1 and the plus sign to 7 2 , 
respectively. 

The integration over momentum p can be easily done after extending the upper limit of integration to 00, so that 
one can use the orthogonality of Laguerre's polynomials 

00 

in 4. n<V 2n 2 
dxe- x x a L^(x)K(x) =T(l + a) { I >' 5 mn with x=-%, (4.8) 

n!a! eis 



and obtain electrical conductivity in terms of the sum over the transitions between neighboring Landau levels 
<t(B, T) = e 2 



„2 v 2 f + v d eBT ' 2 
vpvo 2n 2 T 



X 



v ' {uo 2 + Ml + T 2 )(m 2 + M 2 +1 + L 2 ) - 4i^ 2 A 2 (49) 

„_„ . — -rr ^ + M » + r2 ) 2 - 4^ 2 M 2 ][^ 2 + M 2 +1 + L 2 ) 2 - 4„ 2 M 2 +1 ] ' 



00 - 1 
^7 ^cosb 2 ^' 



The sum over n in Eq. (|4.9(l can be expressed via the digamma function -0 as described in Ref. l2.'"il |. and the final 
expression for the electrical conductivity is 

00 

2 «/jfdh^ ( " ,s,r,A(B)) ' ( «°> 

— OO 

where we introduced the function A 

■4KB,r,A(B)) * * ■ (" 2 + ^ + r^)(eB,»-^(^-^ + r> B 



7T 2 (eS) 2 + (2cjL) 2 [ (w 2 - A 2 -T 2 ) 2 +4w 2 L 2 

^(^ 2 -A 2 + r 2 ) inr ,, ^ a 2 + l 2 -^ 2 -^ J (4 ' n) 

and 

Wr Dn , 

« = — + —. (4.12) 

Another representation of Eq. I|4.11|l can be obtained using the series representation of ip function and writing the 
expression in curly brackets in fractions of l/(r 2 + x 2 ), 

. . . , 1 L 2 f 2 ^ + eBoj(u + A) 4. e Boj(uj - A) 

N (4.13) 
A 2 + M 2 + 2u;M n 2uM n — A 2 — M 2 ' ' 



n=l 



(cc + M„) 2 +L 2 ( W -M„) 2 +T 2 



The representation (|4.13|) is particularly convenient for studying the narrow width limit r«T, \feB when we can 
replace the fractions T/(r 2 + a; 2 ) by ttS(x): 



r r 1 



7T {(eB) 2 +4A 2 L 2 



(f|£ J( „ +a) + (f|)!, ( „_ a) 



O/ m2 1 ( 4 - 14 ) 
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where we kept T 2 in the denominators in order to be able to reproduce a smooth behavior of <r(B) and k(B) in the 
limit B — > 0. We are now in a position to study different asymptotic regimes defined by different relations among the 
parameters T, T, [i, B, and A. 



C. Zero magnetic field 



We begin with the limit of vanishing magnetic field (B = 0). Using the large z asymptote of the ip function 

,i \ i 1 1 1 ^ ( 1 



we obtain 



A(u,B = 0,r,A) 



2tt 2 



2z 12z 2 120z 4 



w 2_ A 2 + r2 A 2 + r 2 -^ 2 

1 H : — ; arctan ; — ; 

2Ur V2 2Ur 



Putting A = also in Eq. JOSJ, we get 



A{u>,B = 0,r, A = 0) - o _. 



1 

27 2 



, w 2 + r 2 w 

1 H arctan — 

r 



(4.15) 



(4.16) 



(4.17) 



Limit T -» 



The limit T — » is significantly simplified by the fact that cosh 2 term in Eq. (|4.1U|I can be replaced by the 5 
function 



cr = e 2 aA(iJ l ,B = 0, T, A), 

where A is given by Eq. I|4.16|l . For the case of zero gap A = using Eq. I|4.17|l we obtain that for fj, = 



<r{n = 0) 
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(4.18) 



(4.19) 



and for \fi\ 3> T 



e 2 a 
4tT T~' 



(4.20) 



Comparing Eq. 14.19|) with Eq. (57) of YN, one can see that our result is twice as large, but contains the same prefactor 
a that reflects the fact that the electrical current operator has the component ~ V£>. On the other hand, the result 
of KC is ~ vf/vd due to the fact that they took the current without the V£> component, and it is twice as large as 
Eq. (|4.19|l because KC integrated over the full Brillouin zone instead of the reduced one. Expression l|4.19|l can also 
be compared with its superconducting counterpart ^SC = (e 2 /tt 2 )vf A>a, where i>a is the dSC gap velocity. As 
one can see the numerical prefactor is exactly the same, but there is no va/vf term in use because the electrical 
current in the dSC state is ~ v^. 

Finally one can also compare our /i ^ expression l|4.20|l with Eq. (13) of KC: 



ct(T< D ) ~ , 

2nv D 70 



(4.21) 



where in the Born limit 70 is related to T(oj) via T(ui) = 70 [(w + /j^/Dq]. Substituting the value of T(0) in Eq. (|4.2(JI) 
one can see that it reduces to Eq. (|4.21|l except for the abovementioncd differences in the factors. 
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2. Limit T < V 

In order to obtain low-temperature corrections to conductivities it is convenient to apply the Sommerfeld expansion 



4Tcosh 2 ^ 



— f(ui)du> 



4Tcosh 2 ^ 



dojfi/j,) 



+ 



oo oo 

/• (^-^)" +1 , 1 /■ 

7 4Tcosh 2 ^ /M+ 2 7 



(w - /z) 



n+2 



4Tcosh 2 ^ 



—diof'iii) + O 



T 



n+3 N 



Here the function / should be nonsingular and not too rapidly varying in the vicinity of uj — fi and the into 
the right hand side of Eq. (|4.22|l are 



duj 

4Tcosh 2 ^ 



du{uj - /i) 2 _ TT 2 T 2 



4Tcosh 2 ^ 



dw{oj - fif _ 7ir 4 T 4 



4Tcosh 2 ^ 
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It 



is easy to derive the next to the leading term in T 2 to the conductivity (|4. for A = 0, 



cr = e 2 a 



ir 2 T 2 

A{lo = fi, o, r, o) + —^—a^(lj = fi, o, r, o) 

o 



where the derivatives of A with respect to u> are 

A! u (w,o,r,o) = 



1 / ! jJ- - 



2-k 2 uj 



arctan — , 

cur r 1 



^(w,0,r,0) = — ^ (- arctan ~ - ' 



tt'^lu' 2 \lu r + r 2 

For [i — case Eq. (|4.24|) gives the T 2 correction to the expression 1|4.19[) : 

7T 2 T 2 



a(n = 0) = 



9 r 2 



T<r. 



(4.22) 



grals on 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



3. Limit T > T 

For r — > from Eq. (|4.16[1 we get 

.Afo,, 5 = 0, r, A) = ^^-^^(c 2 - A 2 ). (4.27) 
For the A = case retaining more terms in the expansion of Eq. (|4.17() , we obtain 

^,B = 0,r,A = 0) = i-(M + 1 L). (4.28) 
Substituting Eq. 14.28|) in the expression (|4.1U|) for a we arrive at 

(J = e2a Sr ln ( 2cosh ^)> T » r - ( 4 - 29 ) 

As one can easily see for [i T Eq. (|4.29|) reduces to Eq. (|4.2U|) reflecting the fact that when /x is the largest parameter 
the value of conductivity is not sensitive to the relation between T and V. 

As mentioned before, the magnetic catalysis phenomenon [2l], 12^, 0] implies that the gap A is generated only in 
the presence of an external field. Nevertheless, one can also study the consequences of the gap opening even for B = 
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case to gain a deeper insight to a more complicated case B,A ^ 0. In particular, considering the A ^> T case, we 
obtain from Eqs. jH2Zt and (|4~TU|) tnat 



167T7T 



duj- 



cosh 2 ^ 



— + (a* -A«) 



At 2 - A 2 ^ 2 T 2 A 2 
|A*| 3AZ 2 |Ai| 
T-p< 2Tln2, 

7 / A 

4T cosh £ exp (-- 



H > A, 

iMl = A, (4.30) 

|m|<a. 



It is clear from Eq. (|4.30|) that the opening of the gap A results in the thermally activated behavior of conductivity 
only for \fj,\ < A. As discussed in Sec. IIVD31 this observation remains valid even in the presence of an external field. 



D. Nonzero magnetic field 

There are not so many cases available for analytical investigation for B ^ and we have to integrate numerically 
Eq. 1(4.10(1 with A given by Eq. ((4.11(1 . Nevertheless in a few cases analytical expressions for the conductivity can be 
obtained and here we begin with considering these results. 



1. Limit T->0 



As one can notice, Eq. (|4.18|) is in fact valid even for nonzero B because only the first of term of the Sommerfeld 
expansion (|4.22|l contributes, i.e., 



a = e 2 aA(n, B, T, A), Va*, T ^ 0. 



(4.31) 



2. Narrow width case 

To study the narrow width limit T — > we use the representation (14.14(1 to arrive at 
__„2 r f {eBf 



AttT 1 2[{eB) 2 +4A 2 r 2 ] 

2(eB) 2 n 



cosh 2 A±l£ 



1 

cosh 2 A^i£ 



J (eJ5) 2 + 4(A 2 + 2eBn)T 2 



1 



1 



cosh 



2 VA 2 +2eB?t+^ 
IT 



cosh 



2 x/A 2 +2eB?t-^ 



2T 



(4.32) 



The fact that the conductivity is proportional to the scattering rate T in the limit T — > means that in contrast to 
the zero-field case (see Sec. IIV CJI) it results from transitions of quasiparticles between neighboring cyclotron orbits. 

Using Eq. I|4.32|l one easily gets the asymptotic expression for conductivity in the strong-field limit, \fe~B > 4T for 
A* = A = 



2 r 

that shows that a becomes field independent in the strong field. 



(4.33) 



3. Numerical calculation of electrical conductivity 



To investigate numerically the behavior of electrical and thermal conductivities, and to make comparison with 
experiment, we need to restore all model parameters, such as K, c, ks, vp, and vp> in Eqs. 1(4.10(1 and 14.11(1 . As 
discussed after Eq. 1(B12(1 the prcfactor a in Eq. ((4.10|) is already fixed and one should only substitute T — * fc^T, 
r — ► fiT, and eB — > (hvpvp)/c)eB. It is convenient to measure all energetic quantities in K, which results in the 
following replacement: 



V2e~B 



hvpVp>2eB 



[K] = 63.9. / — v F [eV ■ A]J B[Tes\&] = 4.206 x 10" 6 . —v F [cm/s] JB [Tesla] , (4.34) 
■>'•■■ y vp 
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where in the first equality vp is given in eV • A and in cm/s in the second one. In particular, for vp = 1.5eV ■ A 
[this roughly agrees with the value of t given after Eq. i|2.21[l ] and vp/vpj — 24 using Eq. I|4.34[l we obtain that 
eB -> 200 • K 2 • B[Tesla]. In what follows we use this estimate to compute all numerical expressions. There is a larger 
uncertainty for the value of T that we could choose for our computations. For example, for clean YBCO monocrystals 
the estimated 30] value of the scattering rate due to impurities To ~ 1 — 2K, so that we use the value T = 2K. 

In Figs. 01 and 0] we show the temperature dependence of the conductivity for three different values of the applied 
field at half- filling (Fig. [21 and slightly away from it (Fig. QJ. For /j, = B = the conductivity increases as the 
temperature T grows. However, when the magnetic field is nonzero, the dependence cr(T) becomes non-monotonic 
and there is first a decrease in a(T) as the temperature increases. For /i, B ^ (see Fig. the conductivity slowly 
increases as the temperature grows and this dependence becomes almost flat as the field increases. This tendency 
agrees with Eqs. I|4.32ll and (|4.33|) . In zero field both figures reveal an '"insulator" (i.e., increasing with temperature) 
behavior in agreement with our analytical result H4.29fl . This behavior is a consequence of using a constant value for 
the scattering rate T in our model. The growth of the conductivity with increasing temperature is directly related 
to the increasing number of thermally excited quasiparticles. The choice of temperature-independent T might be 
reasonable in the narrow low-temperature region. It is essential, however, that the observed decrease and flattening 
<r(T, B) in the nonzero external field are due to the assumption T(B) — const discussed after Eq. i|3.8fl and are not 
related to the r(T) = const approximation just mentioned. 

Comparing Figs. and 01 one can see that in zero field the increase of \fi\ (opening of the pockets on the Fermi 
surface) leads to the increase of the conductivity. On the contrary, in the presence of the external field, the conductivity 
decreases as the pockets on the Fermi surface open. One can notice that the value of conductivity a at fi = T = 
is field independent. This is due to the fact that for oj = 0, the function A [see Eq. (14. 1 1|> and Fig. ^] becomes field 
independent. 

Using Eqs. I|4.10|l and l|4.11|l it is also possible to investigate the case of A ^ 0. The main results for B ^ can 
be foreseen from Eq. (|4.30(l and are the following. When T, \fj,\ < A the behavior of a becomes thermally activated, 
i.e., governed by the factor ~ exp(— A/T). In contrast, for A < \/j,\ the behavior of cr(T) is rather similar to the case 
of A = 0. This reflects the fact that the gap A induced by the magnetic catalysis, in contrast to a superconducting 
gap, is not tied to the Fermi surface and there are gapless excitations for > A. 

In Figs. [5] and [S] we show, respectively, the dependence of <j(B) for three different values of T and fixed fi and 
for three different values of fi at fixed T. Both figures show that a decreases as a function of B. This decrease 
reflects the fact mentioned before Eq. I|4.9|l that in the presence of magnetic field only transitions between neighboring 
Landau levels contribute into electrical (and thermal) conductivity. A further increase of the field and entering in the 
strong field regime, v 'eB > AT suppresses the transition between Landau levels, so that as one can see from Fig. [S] 
the conductivity o~(B) becomes field independent. Furthermore for T = 10K the value of normalized conductivity is 

0.15, in agreement with Eq. (|4.33|1 . which is derived under the assumption T -C T and taking into account only 
the transition between the lowest levels. In Fig. Elwe observe Shubnikov-de Haas oscillations of the conductivity that 
are due to the Landau-level crossing of the Fermi pockets (an analytical treatment of the magnetic oscillations in the 
conductivities will be presented elsewhere). 

V. THERMAL CONDUCTIVITY 
A. General expression for thermal conductivity 

Thermal conductivity can be calculated from the the energy current-current correlation function 



U E Jm =-j dre^<T T jf (r,0)jf (0,0)), 


and the correlation function of energy current with the electrical current 



= -J dre^(T T j^(r,0)MQ,0)) 


using a thermal Kubo formula [27|, US] 

_ _ 1 imnjP(» + »o) 1 imng g (n + iO) _ 1 [imnjf(» + »o)] 2 



(5.1) 



(5.2) 



(5.3) 
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Here Jl EE (fl + iO) = U EE (iQ, -> Q, + iO) is the longitudinal polarization [see Eq. (|B14|) ] and S(CL,B,T) is the 
thermopower 



S(Sl) 



1 lntfl| c (0) 

ri m ng c (fi)' 



(5.4) 



with the longitudinal nf c (O-H0) = n BC (iO -> tt + iO). The term with the thermal power S ensures that the energy 
current is evaluated under the condition of vanishing electrical current (see also Ref. [2?! )■ Usually, for T <C [i this 
term is considered to be unimportant because it is ~ T 2 / fj, 2 times less than the first term of Eq. (|5.3|) . It is also zero 
for /i — 0, but nevertheless its contribution is important for the case of interest, < T. 

The structure of the energy current (|2.18|) is similar to that of the electrical current l|2.14[l . Hence, the calculation 
of the corresponding polarization functions (|5.1|l and l|5.2|) is almost identical to calculating the bubble l|4.1|l . Using 
the general result (|B16|I to the case of interest: g = g' = to + 0/2 — fi, we arrive at 



Imn^Q) _ 7T 



d 2 P 



dio- 



tanh "~£+ n - tanh ^ (to - \i + ft/2 



T 2 Q v F v D J (2tt) 2 J fl \ T 

— OO 

x [4tr[A z ,(w,p) 7 M I 3(w + fi,p)7 1 ] + v 2 D tr[A D (uj,p)j 2 A D (uj + n,p)j 2 ]] , 



(5.5) 



and for g = to + 0/2 — /i, g' = e 
lmnf c (0) _ en 



Tfl 



V F VD 



d 2 p 



tanh^±^ 
dui — — 



2T 



to - n+ £1/2 



T 



(5.6) 



x [v 2 F tv [A d {uj, p^Aoiuj + Q, p) 7 1 ] + v 2 D ti [A D (u, p)j 2 A D {u + ft, p) 7 2 ]] 



Similarly to the case of the electrical current, one could also obtain the energy current operator from the nodal 
Lagrangian i|2.38[) in the external magnetic field as done in Ref. [23j| . 



■ E lV F i - 1 Q a - 1 

Jx = —KXsl OtXs ~ OtXsl Xs 



. E _ lv D 



(5.7) 



to derive the expression l|5.6(l . The current l|5.7|) corresponds to the current used by YN. 

The calculation of thermal conductivity from Eq. 1)5. 6|) follows exactly the same route as for the electrical conduc- 
tivity in Scc. lIVBl and finally we arrive at 



K (B,T) 




4Tcosh 2 ^ 



—A(lo,B,T,A(B)) 



2 

e a 



a(B,T) 



LU — \l 



T 4Tcosh 2 ^ 



(5.8) 



where A is the same function 14.11)1 as for the electrical conductivity and a is given by Eq. (|4.12)l . Thus we are ready 
to study the thermal conductivity as a function of T, T, fi, B, and A. 



B. Zero magnetic field 



To study the B = case we substitute the expression for A(B = 0) given by Eq. (|4.1(i|) into the general expression 
for thermal conductivity (|5j 



Limit T -» 



Corresponding to Eq. (|4.18(l limit T — > is given by 



K 7T 

T ~ ~3 



olA(jjl,B = 0,T, A). 



(5.9) 
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Then using Eq. 1(4.17(1 for // = A = 0, we obtain 



k(h = 0) 



(5.10) 



T 3 

The factor a in Eq. (|5.1UI) is present also in the DDW 0, @ and dSC 01 cases where the same expression for 
the thermal current has been used. The overall numerical factor 1/3 in Eq. 1(5.10(1 is the same as in Ref. ^3 an d 
it is different from Refs. 0, 0] in the same way as the corresponding factor in expression ((4.19(1 for the electrical 
conductivity. 



2. Limit T < V 



Using the Sommerfeld expansion ((4.221) we can also derive next to the leading term in T 2 to the thermal conductivity 
for A = 



k a 
T ~ 3 



(5.11) 



where the derivatives of A are given by Eq. (|4.25|) . For fi = 0, since there is no contribution from the thermal power 
term II , Eq. (|5.11|) reduces to the expression 



k(/j, = 0) a 
f = 3 



15" r 7 



t< r. 



(5.12) 



3. Limit T > T 

For T ^ and A = substituting Eq. (|4^H| into the expression H5.8(l for k, we obtain (see Appendix EJ) 

" 7T [J, 7T 3 T 2 ' 

12T ~ 36 I> 



K 



(5.13) 



where the first term ~ [i of K arises from the first term of Eq. I|5.8|l and the second term ~ l//i originates from the 
thermal power term II^. Thus, as one would expect, the thermal power contribution into the thermal conductivity 
is important only for \fj,\ < T. 

For the case /j, = A = using Eq. 14.28(1 we arrive at the expression 



n(p = 0) 
T 



9C(3)T ln2T 

4 r + ~Yt 



r«r 



(5.14) 



where (,(z) is the Riemann zeta function. Finally using the representation ((4.27(1 . which is valid for // = we can get 
the thermal conductivity for nonzero A |23|: 



k(/x = 0) 



a 



87rT 3 r 



A 



kl(^ 2 -A 2 ) 
cosh 2 ^ 



A 2 

( 

tiTT 



exp 



A>T. 



(5.15) 



C. Nonzero magnetic field 

As for the electrical conductivity there are not so many cases available for analytical investigation for B / and 
we have to integrate numerically Eq. 1(5.8(1 with A given by Eq. ((4.11(1 . There are still a few cases when the analytical 
treatment is possible and we begin by looking at them. 



1. LimitT^Q 



As one can check, Eq. (|5.9|) is valid even for nonzero B, so that 

2 

l = ^-aA( f M,B = 0,T,A), V/*, r ^ 0. 



(5.16) 
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2. Narrow width case 



Using the narrow width representation l|4.14|l for A, we arrive at the following result, 



K 
f 



(eB) 2 



4ttT 3 2[(eB) 2 



E 



(A + M ) 2 , (A - p) 2 



• 4A 2 r 2 ] 

2{eB) 2 n 



cosh 2 ^ 



cosh 2 ^jf^ 



(eB) 2 + 4(A 2 + 2eBn)T 2 



(VA 2 + 2eBn + (V A 2 + 2eBn - n) 2 



cosh 



2 VA' 2 +2eBn+Li 



2 ^A' 2 +2eBn-L 



cosh 



T ' 



2T 2T 

where ki is a term originating from the condition of vanishing electrical current. As was discussed above, 
for [i = and can be neglected for 3> T. In particular, for /i = A = in the strong field limit, \fe~B 
obtain from Eq. (|5.17|l the expression 



K 
f 



8TeB 



(5.17) 



it is zero 
> 4T we 



(5.18) 



3. Numerical calculation of thermal conductivity 

In Figs. [7] and [H] we present the results for the thermal conductivity that were obtained for the same values of the 
model parameters as the data for electrical conductivity shown in Figs.|3and0J Comparing these figures, one can see 
that the behavior of the thermal conductivity is rather similar to the behavior of electrical conductivity and to see 
more subtle differences we have to consider the temperature dependence of the Lorenz number as done in Sec. IVII 

It is also interesting to investigate numerically the contribution from the second term of Eq. (|5.8(l to the thermal 
conductivity. This is done in Fig. EI where one can clearly see that when the temperature grows and T becomes < \fi\ 
the second term of Eq. 15.8fl makes an important negative contribution to the thermal conductivity. 

In Figs. HOIandllll we present the results for the dependence of thermal conductivity on the magnetic field B that 
were obtained for the same values of the model parameters as the data for electrical conductivity shown in Figs. [S] 
and[|)l 



VI. WIEDEMANN-FRANZ LAW 



After deriving in Sees. II Vl and IVl the expressions for electrical and thermal conductivities in the various limits we 
are ready to consider their implications for the WF law. 



A. Limit T — > 0, arbitrary field, and chemical potential 

For finite T, arbitrary /i, B, and A in the limit T — > one can see from Eq. (14.3111 [see also Eq. I|4.18|l for B = 0] 
and H5.16H [see also Eq. (|5.9() for B = 0] that WF law is maintained: 



K 

of 



TT^k 

Ye 



(6.1) 



This result can be understood from the qualitative arguments given, for example, in Ref. The expressions for the 
electrical, Eq. Q4. !()[! . and thermal, Eq. Ij5.8|l . conductivities calculated within the bare bubble approximation are very 
similar in the limit T — * when the second term of Eq. (|5.8J) vanishes. Both expressions contain the function A(u>) 



multiplied by the derivative of the Fermi distribution f(u>) 



n(cj — /i), with the only difference that the thermal 



conductivity is also multiplied by the factor ~ (lu — fi) 2 . These prefactors f{u) and g(uS) = —(lu — (i) n' F (uj — fx)/T 2 
are shown in Figs. Upland IT31 When the temperature T goes to zero, the 5-like spikes of f(u>) and g(v) occur at the 
same value u> = fi and the WF law is maintained. 



B. Zero magnetic field at T / 

For (j, = and T » T considering Eqs. (|5.14|) and 14.29fl we obtain that the WF law is violated, 

k 9C(3) k% 
L -^f-2^2^- 2 - 37L °' 



(6.2) 
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while for p 3> T Eqs. I|5.13|) and (|4.29(l show that the Lorenz number has its usual value 

L = L . (6.3) 

These two results can also be understood qualitatively using Figs. El an d El For \fi\ ^> T > the electrical 
conductivity is determined by energies u> w p, = 0, while the thermal conductivity is determined by the energies near 
u> sa p ± /csT. If the function ^4(w) does not vary appreciably over the energy range p — fc^T to p + fc^T, the Lorenz 
number may still be ~ 1 even for finite T,T <C //, as follows from Eq. I|6.3() . If /i becomes of the same order of T, the 
dimensionless values p/T and p/T ± 1 are apart and the Lorenz number deviates from 1 as seen from Eq. 1)6.2(1 . 
For nonzero A > T and (i = 0we get from Eqs. I|4.3l)[) and l|5.15l) that 

k 2 A 2 

^=?^ (6-4) 

i.e., L > Lo for A > (n/y/S)T. This dominance of the thermal conductivity can be easily understood from the fact 
that when the gap opens, the electrical conductivity diminishes more strongly because it is determined by the energies 
lj sa p near the Fermi surface. 



C. Nonzero magnetic field at T / 

For p = A = in the strong field limit, \fe~B > AT from Eqs. (|4.33l) and l|5.18|l we obtain that the Lorenz number 
becomes field dependent: 

32eB _ vSE 

L = -j^-e t . (6.5) 

In the presence of magnetic field the function A(lo) (see Fig. I14|) varies appreciably over the energy range u> — ksT 
to uj + ksT, so that it is nonzero for u> = and very small at u> = ±ksT. This feature of A(uj), along with the small 
values of p, can produce rather strong violation of the WF law. 

Our qualitative arguments are indeed confirmed by Figs. IT31 and 1161 where we show the temperature dependence of 
L(T) at half-filling (Fig. I15fl and away from it fFig. 1161) . Both these figures are computed on the basis of Eqs. (|4. 1 U|) . 
(|5.8|) . and l|4.11|l . For p = the value of L always goes down from Lq as the temperature increases, then L increases 
crossing Lo, and finally it goes to its zero-field value l|6.2() . One can see from Fig. El that f° r nonzero /i the deviations 
of L(T) from L become less pronounced (see also the discussion of Fig.lT7l). 

It is also instructive to investigate the importance of the second term of thermal conductivity, Eq. I|5.8|) . for the 
Lorenz number calculated with and without this term (as already done in Fig. for the thermal conductivity itself). 
These results are given in Fig. 1171 and they show that the second term of thermal conductivity is indeed crucial to get 
the correct answer. Note that for the case of nonzero /z the situation with the WF law is more complicated because 
the value of /x may coincide with 10, corresponding to the maximum of A(lu). In this case L(T) will firstly increase, 
as one can see in Fig. E|f° r B = 12T. 



VII. CONCLUSIONS 

Let us compare the results derived in the present paper with the experimental results obtained in Ref. It is 
still unknown whether the DDW state exists and/or plays an important role in the electron doped compound used 
in Ref. QjjSO that making this comparison would help to address these questions. First, we observe that while the 
results of [j| suggest that the WF law is violated at T — > 0, there is no violation of the WF law in this limit in the 
DDW scenario of the pseudogap. Since in Fig. 3 of Ref. Q the electrical conductivity is a constant, the line n e (T)/T 
directly represents the normalized Lorenz number L(T)/Lq. For finite temperatures there is then some similarity 
between Fig. 1151 for Fig.EJf° r M 7^ case) and Fig. 3 of Ref. where as T increases the thermal conductivity crosses 
from the region with K e /T < Lq/pq to the region with K e (T)/T > Lq/po resembling the character of the WF law 
violation seen in the experiment. 

However, in the experiment, the electrical conductivity is flat while the thermal conductivity changes significantly in 
a subkelvin range. On the other hand, theoretical calculations show that the electrical (see Figs. |3 and 0} and thermal 
(see Figs. and |5J) conductivities vary simultaneously in a range of 10 K. Note that the wider range of temperatures 
can be probably attributed to the fact that the values of the model parameters we took for the numerical estimates 
are more appropriate for the hole-doped compounds. 
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It is obvious from Fig. 1151 that such a behavior of L(T) is due to the presence of the magnetic field. This confirms 
our claim that to interpret theoretically the experiment in Ref. pj one should take into account the influence of the 
external field. 

Our main results can be summarized as follows. 

(1) We have obtained analytical expressions for electrical conductivity l|4.10[l and thermal conductivity (|5.8() in the 
DDW state in presence of an external magnetic field. 

(2) We have established a correspondence between the expression for the electrical conductivity l|4.3[l in the DDW 
state written in terms of the generalized velocity V(k), Eq. I|2.15|l . and the dc conductivity derived by Langer |28|] in 
the bubble approximation with the vertex l|4.4|l . 

(3) We have shown that in the DDW system in the presence of impurities the WF law holds in T — > limit for an 
arbitrary field B and chemical potential \x. This is checked within the bubble approximation, i.e., not including the 
impurity vertex. The influence of the impurity vertex on the properties of the DDW state in zero field at half-filling 
can be considered using the heuristical arguments of Durst and Lee 01 f° r dSC state. Since in the dSC state the 
thermal (and spin) currents are proportional to the group velocity, they can relax through either intranode scattering 
or scattering between nodes. As a result, the different types of scattering play nearly the same role and therefore 
vertex corrections do not contribute to the thermal (and spin) conductivity. This is not the case of the electrical 
current in the dSC state, because it depends only on the Fermi velocity. Thus, the electrical current can relax more 
effectively via scattering from node to node than it can via scattering within a single node. This difference is taken 
into account by considering vertex corrections which modify the bare bubble expression for the electrical conductivity. 
Since the electrical current in the DDW state is also proportional to the group velocity, the corresponding impurity 
vertex should not contribute to the electrical conductivity. Moreover, the thermal current in the DDW state is also 
proportional to the group velocity, so that one can rely on the same arguments about the impurity vertex as for the 
dSC case. Thus we do not expect that in the DDW state the impurity vertex can be a source of the WF violation 
at T — > 0. These observations should be confirmed by detailed calculations of the impurity vertex corrections, which, 
however, are beyond the scope of the present paper. 

(4) For finite temperatures T < the WF law violation is possible and in zero field the thermal conductivity 
dominates over the electrical conductivity, i.e., L(T)/Lq > 1. 

(5) For T < \fi\ in the nonzero field the WF law violation becomes even stronger than in zero field and depending on 
the temperature both regimes L(T)/Lq <C 1 and L(T)/Lq > 1 are possible. 

(6) For T <C there is no WF violation even in the presence of magnetic field. 

Finally, we would like to stress that the results of the present paper should be applicable not only to the DDW 
state, but to a wider class of theories. Here we started from the DDW Hamiltonian (|2.1|l and, to simplify the problem, 
approximated it by the QED3 Lagrangian l|2.38[) . There is, however, a number of systems, e.g., pyrolitic graphite |24| . 
that can also be described by the Lagrangian (|2.38|) . so that our results would also be relevant for them. 
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APPENDIX A: NODAL FERMION GREEN'S FUNCTION IN AN EXTERNAL MAGNETIC FIELD 



In the coordinate space, the fermion propagator has the following form in the proper-time representation 

-1 



S(x,y) = (W + A) x {x 



A 2 + D 2 



y) = -i(iD + A) x / ds(x\ exp[-is(A* + D z )]\y), 



(Al) 



where D — 7" ' D u and the "long" derivative D v is given by Eq. 12.37(1 . The matrix element (x\ exp[— is(A 2 + D 2 )]\y) 
can be calculated using the Schwinger (proper time) approach 22] (see also a pedagogical overview in Ref. [3 ) . The 
main idea of this method is based on the interpretation of (x\ exp[— isD 2 ]\y) as the coordinate representation of the 
proper-time "evolution operator" U(s) = exp[— iHs], where we have introduced the "Hamiltonian" 



H = D 2 = D 2 



-<j pv F 
2 p u > 



F pu = d p A v - d v A p , a 



pv 



(A2) 
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The matrix element of the evolution operator can be evaluated Jl9|. for example, using either operator or functional 
integral formalism: 



(x\U(s)\y) 



8(7TS) 3 / ; 



■ exp 



u j dt v A»(0 + % -{x - y)f(s)(x -y)- L(s) + i^a¥s 



with 



f(s) = eFcoth(eFs), L(s) = -trln 



/ sinh e¥s 
\ e¥s 



(A3) 



(A4) 



where we have used matrix notations, e.g., F pv = (F) p „ and the integral with A v is calculated along the straight line. 
This operator differs from the four-dimensional version of (x\U(s)\y) |19|| in the power dependence of the proper time 
and in the numerical prefactors. Plugging (| A3|) into IjAljl one obtains 



5(x,y)=exp lie / d&A"® 



x exp 



oo 

e - l37 r/4 r ds 

J ^ 




8tt 3 / 2 



^--Y[f{s)+e¥] pv {x-yY 



(A5) 



-iA 2 s + -(x - y)f{s){x - y) - L(s) + i-aFs 



Finally, to perform calculations for the case of interest one should evaluate f(s), L(s) and exp[i(e/2)erFs] for F 
containing purely magnetic background field B. The corresponding expressions for these values are 



f(s)pu = - 

exp[-L(s)] 
exp 



9 P , 



„ + ^f(l-e5 S cot(eB S )) 
B z 



eBs 



(A6) 



i— a¥s 
2 



sineSs 
= cos eBs + j 1 sin eBs. 



In the Matsubara frequency-momentum representation this leads to Eq. (|3.2|l . 

It is convenient to write down an alterative form of the fermion propagator in a magnetic field as a sum over the 
Landau-level poles. In deriving it we follow Refs. |2lJ,|25|. Introducing the shorthand notations we rewrite Eq. i|3.2|) 



as follows, 



ds exp 



s [a + p 



, tanh(e-Bs) 
eBs 



x [(6 + ctanh(| eJ Bs|)) (1 + titanh(|eBs|))] 



where 



a = A 2 - (iui) 2 , b = iujj° ~ p^ 1 - p 2 J 2 + A, 

c = —i(p2j 1 — pi7 2 )sgn(e_B), d = — i7 1 7 2 sgn(eB). 

Using the identity tanhx = 1 — 2cxp(— 2a;)/[l + cxp(— 2a;)], the relation [32| 

(l-zr^exp )=j2K(x> n , 

^ ' n=0 

2 

with z = — exp(— 2|e£?|s), x = 2y^gj, we arrive at the expression 
S^ito, p) = — exp 



z < 1, 



(A7) 



(A8) 



(A9) 



x ds exp 



\eB\ 
2p 2 z 



eB z- 1 



{b + c) (1 + d) + (4cd + 2c + 2bd) -^—^ + 4cd - ^ 2 



(A10) 
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Then using the identity L" 1 (x) = L™(%) — and the definitions L n = lP n , L" 1 =0we obtain 



S(iu>, p) = — exp 



eB 



ds exp(— sa) 



(b + c)(l + d)J2L n {x)z n 



(- (b + c)(l + d) + 2c + 2bd) J2 L n -x{x)z n + 4cdJ2 2£-i(a0* 



n=0 



n=0 



exp 



eB 



J dsexp(-sa) < (A + iw 7 °) 



(1 - sgn( e S) l7 1 7 2 ) L n {x)z n - (1 + sgn(eB)i 7 V) £] L n -x{x)z n 



n=0 



n=0 



+4(p 1 7 1 +P27 2 )E L n-i(^ n ' 

n=0 ) 

Finally, integrating over s we arrive at Eq. 



(All) 



APPENDIX B: GENERALIZED POLARIZATION BUBBLE AND NODAL APPROXIMATION 

The calculations of electrical, thermal, and even spin conductivities are quite similar for both DDW and dSC 01 
nodal systems. Thus instead of repeating the same calculation several times and to underline the similarities and 
differences between DDW and dSC cases, it is rather convenient to define a generalized polarization tensor H" that 
depends on the generalized coupling parameters g, g' 



9,9 



iu> + iO/2 
u + 0/2 - fi 



(Bl) 



Two lines in the definition of g, g' correspond to the Matsubara and real frequencies, respectively, and the origin of 
the chemical potential /1 in the second line will become clear later. The generalized bubble is 



n^(»fi) = 2T 



d 2 k 



RBZ 



(27T) 



2 E 



fla'tr[G(zw,k)V r Q (k)G(iw + in,k)V^(k)], 



(B2) 



where G(£cJ,k) is the Green's function l|2.2U[l . but evaluated in the external field, and V(k) is generalized velocity 
given by Eq. (|2.15l) . The integral is over the reduced Brillouin zone (RBZ) and the factor 2 before the integral is due 
to the spin degree of freedom s. To compare this expression with dSC case [see Eq. (Al) of Ref. 0] we n °te that 
the summation over the spin degree of freedom is already included in the Nambu formalism. 
Using the spectral representation for the fermion Green's function 



re x w ? ^i> k ) A 
Gucj + /i,k) = / ajJJx, 

J IUJ + fl — Ml 



where A{uj\, k) is the spectral density, we arrive at 

d 2 k 



nM) = 2 



RBZ 



{2-Kf 



dut / du 2 tT[A(u Ju k)V a (k)A(Lu 2l k)V p {k)}R 



with 



R = TY gg'-. , l — . 



(B3) 



(B4) 



(B5) 



Since the intermediate results differ depending on the frequency dependence of the coupling parameters g, g', 
we consider the frequency-independent and frequency-dependent couplings separately. For g = g' = e (frequency 
independent coupling), evaluating the sum and then continuing iQ — ► 17 + iO we get 



2 nf(wi — fi) — np{u2 — fi) g 2 tanh U2 2T ^ — tanh "^t^ 

R = g 



u>i — U02 + n + io 



u>i — uj 2 + + iO 



(B6) 
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For the frequency-dependent coupling g = g' = ito + ifl/2 the evaluation of the Matsubara sum i|B5(l gives [3£ 
(iuj + ifl/2) 2 (u)% — [i + iQ/2) 2 riF(ui — fi) — {0J2 — M — ifl/2) 2 np{^2 — m) 



(iuj + fi — u>i)(iu! + fi + iQ — L02) 
Then, continuing to real frequencies iQ — > f2 + iO we obtain 



UJl — Ll>2 + i£l 



ImR= -(uj 1 - fi + n/2) 2 



tanh 



u>i — fi u)\ — fi + VL 



2T 



tanh ■ 



2T 



8{loi — uj 2 + fi). 



(B7) 



(B8) 



Comparing the imaginary part of Eq. (|B6|I with Eq. I|B8|) one can see that the latter equation can be obtained from 
the former by the direct substitution of g = g' = uj + fi/2 — fi. This explains the difference between the upper and 
lower lines in Eq. (|B1(I . 

Similarly, evaluating the Matsubara sum (|B5|I with g = e, g' = iu> + Klj2 we get 



R = T y^ e(iuj + itt/2) 

4-^ {iui + fi — LUi)(iuj + ji + iil — 0J2) 



e(wi — fi + ifi/^^F^i — fi) — e(uji — fi — i£l/2)np(uj2 — /-^) 

UJl — Ll>2 + i£l 



so that 



ImR = -e(uj 1 - fi+ fi/2) 



tanh 



2T 



tanh 



uji — /i + f2 



2T 



- w 2 + fi). 



(B9) 



(BIO) 



Finally, we obtain for the imaginary part of the tensor polarization 11^ (Q + iO) the following expression 

oo 

d 2 k 



Imn^(fi + iO) =tt 



RBZ 



(2tt) 2 



dujgg' 



to — ft uj — ft + Q 

tanh — — tanh 



2T 



2T 



tr [A{u>, k)V a (k)A(w + fi, k)V/j(k)] . (Bll) 



Looking at Eqs. (|B2(I and IjBlljl one can notice that the Green's function G and the associated spectral density 
A in the external field are in fact unknown. Instead of considering these functions, we have constructed in Sec. IIIII 
the Green's function (|3.3I) for the linearized nodal Lagrangian l|2.38[l that is valid in the vicinity of the four nodal 
points. Thus, in Eq. (|B2() we replace the integration over the reduced Brillouin zone by the integral over the k-space 
surrounding each node and sum over the four nodal subzones: 



RBZ 



d 2 k 



2 ^ 

3=1 



(27T) 2 



2 ^ 

3=1 



d 2 p 



(2it) 2 vfvd Akvfvd ~~ 



E 



PO 27T 

pdp 



3=1 o 



de 

2t? 



(B12) 



where pi = vpk x = pcos#, p 2 = vpky — psmQ, p = \/ p\ + pf) = ^/t^VfVd / (2a) and the local nodal coordinate 
systems (k x ,k y ) are shown in Fig. [3 Note again that, comparing with fljjj], an extra factor 1/2 appears due to the 
fact that the original integral is over the reduced Brillouin zone. 

The advantage of the scaled variables pi and p2 is that they correspond to the "relativistic" case vf = Vjj — c = 1. 
Moreover, we can use the spectral function Ad in Eq. (|3.(i[) inside the integral over p, so that all necessary factors 
with vf and vd are already outside the integral i|B12() . We should only provide the rule that allows us to restore the 
model parameters coming along with the magnetic field in the final result |34| : eB — ► (HvfVd /c)eB. 

We replace the generalized velocities V(k) in Eq. (|B11|) by their values Vk=N on the Fermi surface at half- filling 
(fi = 0) . Note that for the DDW case this approximation is more severe than for the cZSC case because for fi ^ the 
true Fermi surface does not coincide with the Fermi surface at half-filling. This approximation puts some restrictions 
on the values of fi, so that we cannot move far away from half-filling. The bubble (IB lip takes the form 



ImlT^ (fl+iO) = 



2vfVd 



E 

3=1 



d 2 p 



duigg' 



tanh ■ 



2T 



uj — \i + Q 
2T 



tanh - '- tr [A D (w, p)j°V a A d {uj + fi, p)7%] , 

(B13) 

where the spectral density l|3.6[) includes the broadening of the spectral lines (|3.8|l due to impurities and we inserted 
the 7° matrix [see Eq. 1)3.9(1 ] . Note that in Eq. (|B3|) we used only the translationary invariant part (|3.2|l of the Green's 
function (|3.1|) calculated in the external field, since its translation non-invariant part cancels out when substituted in 
the bubble II. 
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Since we are interested only in the longitudinal conductivities and because the system is isotropic, we can define 
the longitudinal polarization function IT" as follows (the sum over dummy index is implied): 

n^(n) = in^(fi). (bh) 



Then evaluating the sum over nodes in Eq. (|B13|I by using the identities 

4 4 

J2 v l* v l0 = 2v l S <*0> J2 v li )v i'f> = 2v M' e * e <*P ('^0, v, = {v F ,v D }, (B15) 

3=1 3=1 



with ejj being antisymmetric tensor, we arrive at the final result 



oo 

ImIP*'(ft + iO) = — I I duj 99' 

vfvd J (2ny J 



u> — n lo — fi + n 

tanh tanh 

2T 2T 



— OO 

x (w|tr[A I? (a;,p)7 1 A D (w + fi,p)7 1 ] + v 2 D tv [A D (u, p)-f 2 A D (uj + Q, p) 7 2 ] ) . 



(B16) 



Since both the electrical current (|2.14|) and thermal current l|2.18|l have two terms, a "Fermi" term proportional to 
vp and 03 and a "gap" term proportional to and er 2 , our calculation of the current-current polarization function 
(|B2|I is in fact similar to the calculation of the thermal current bubble for the dSC case done in Ref. ^Tfl • Evaluating 
the polarization function (|B13|) we obtain four bubbles: Fermi-Fermi, Fermi-gap, gap-Fermi, and gap-gap. However, 
since vj?lvj) at each of the gap nodes, the two cross terms cancel, as reflected in the second identity in Eq. 1|B15(I . 
Therefore both electrical and thermal longitudinal conductivities have two terms: the Fermi term with the velocity 
Vp and the gap term with the velocity v 2 D . Thus, finally we approximated the polarization bubble for the Hamiltonian 
(12. 1|) by the polarization bubble for the nodal Lagrangian l|2.38[) . 

APPENDIX C: SOME INTEGRALS 

Here we derive the expression (|5.13|l for k. It follows from the general expression (|5.8|) with A given by Eq. I|4.28|l . 

The corresponding integrals are easily evaluated in terms of polylogarithmic functions, Li„(2;), |35j: 



J = 



dujluj — [i) 2 \u)\ 



—00 



cosh 2 ((w-/x)/2T) (CI) 
= 8TV ln(2cosh( A1 /2T)) + 16T S fi[Li 2 [-e^ T ] - Li 2 [- e -^ T ]] - 24T 4 [Li 3 [-e^/ T ] + Li 3 [-e^ /T ]] 

and 

00 

1 = I Sr^wL = 8T 3 [Li 2 [-e-^] - Li 2 [-e^]] - 87^ ln(2 cosh( M /2T)). (C 2) 
J cosh (w — u) 21 



- DC 



These rather complicated expressions for / and J can be expressed via elementary functions in the limit /i>T. We 
begin with the more simple case of I. Using the identity [^3 

Li 2 (z) = -Li 2 (l/z)-iln 2 (-z)-y, z£(0,l), Li 2 (0) = 0, (C3) 
we obtain that in the limit /i T 

Li 2 [-e"^] - Li 2 [-e^] = 2Li 2 (-e^) + \\n 2 {e^ T ) + £ * + (C4) 
Thus the final expression for I reads as 

4^ 2 



I - —T d , /i » T. (C5) 
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To simplify J, in addition to Eq. (|C3|) . we use the corresponding identity for Li^^z) [35l| : 

Li 3 (z)=Li 3 (l/z)-iln 3 (-z)-^ln(-z), (0,1), Li 3 (0) = 0, (C6) 

6 6 

and for fi>T we obtain that 

j ~ ^-T 3 », M » T . (C7) 
Putting Eqs. l|C7j) and l|C5j) along with Eq. (|4.2flj) into Eq. H5.8[) . we arrive at the final result l|5.1H|l . 
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FIG. 1: Schematic phase diagram of cuprates within the DDW scenario. 
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FIG. 2: Fermi surface in the DDW state. The hole pockets are centered around the points (±7r/2, ±7r/2). We choose the local 
nodal coordinate systems as shown. 
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FIG. 3: The normalized conductivity a/a(fi = T = B = 0) as function of temperature, T, for three different values of magnetic 
field B at half-filling, fj, = 0. The value of a(/i — T — B — 0) is given by Eq. UTTUl l. 
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FIG. 4: The normalized conductivity a/a(fi = T = B = 0) as function of temperature, T, for three different values of magnetic 
field B away from half-filling, [i = 20K. 
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FIG. 5: The normalized conductivity a/a(fj, = T = B — 0) as function of field, B, for three different values of temperature T 
at half-filling, fi — 0. 
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FIG. 6: The normalized conductivity o"/<r(/x = T = B — 0) as function of field, B, for three different values of chemical potential 
(j, at T = 2K. 
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FIG. 7: The normalized thermal conductivity k/k(h, B = 0,T — > 0) as function of temperature, T, for three different values of 
magnetic field B at half-filling, /i = 0. The value of B = 0, T — » 0) is given by Eq. 15.101 , 
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FIG. 8: The normalized thermal conductivity k/kSji, B — 0, T — » 0) as function of temperature, T, for three different values of 
magnetic field B away from half-filling, /i = 20K. 
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FIG. 9: The normalized thermal conductivity k/k((J,, B — 0,T — » 0) as function of temperature, T, for three different values of 
magnetic field B away from half- filling, n = 100K. The dotted lines are calculated without the second term of Eq. 15.81 that 
originates from the condition of absence of the electrical current in the system. 
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FIG. 10: The normalized thermal conductivity n/ft(fi,B = 0,T — > 0) as function of field, B, for three different values of 
temperature T at half-filling, fi = 0. 
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FIG. 11: The normalized thermal conductivity k/k((j,,B = 0, T — > 0) as function of field, B, for three different values of 
chemical potential \i. 




FIG. 12: The expression for electrical conductivity contains the factor f(u>) = —n' F (uj — /i) (the case n — is plotted) and thus 
samples A(u>) in the immediate vicinity of the Fermi surface. 
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FIG. 13: The expression for thermal conductivity contains the factor g(to) = — n' F (uj — — n) 2 /T 2 (the case \jl = is plotted), 
and thus measures A(lo) in the immediately below and above the Fermi energy. 
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FIG. 14: If A(u>) varies appreciably over the energy range from n — ksT to \x + ksT then the Lorenz number deviates from its 
usual value Lq. 
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FIG. 15: The normalized Lorenz number L/Lo as function of temperature, T, for three different values of magnetic field B at 
half-filling, /i = 0. 
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FIG. 16: The normalized Lorenz number L/Lq as function of temperature, T, for three different values of magnetic field B 
away from half- filling, /i = 100K. 




FIG. 17: The normalized Lorenz number L/Lo as function of temperature, T, for three different values of magnetic field B 
away from half-filling, fj, = 50K. The dotted lines are calculated without the second term of Eq. 15.81 that originates from the 
condition of absence of the electrical current in the system. 



